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Abstract. We report on a sensitive search for mid-infrared molecular hydrogen emission from protoplanetary disks. We 
observed the Herbig Ae/Be stars UX Ori, HD 34282, HD 100453, HD 101412, HD 104237 and HD 142666, and the T Tauri 
star HD 319139, and searched for H2 - 5 (2) (7 = 4 - 2) emission at 12.278 micron and H 2 - S (1) (7 = 3 - 1) emission 
at 17.035 micron with VISIR, ESO-VUT's high-resolution mid-infrared spectrograph. None of the sources present evidence for 
molecular hydrogen emission at the wavelengths observed. Stringent 3cr upper limits to the integrated line fluxes and the mass 
of optically thin warm gas {T = 150, 300 and 1000 K) in the disks are derived. The disks contain less than a few tenths of 
Jupiter mass of optically thin H2 gas at 150 K, and less than a few Earth masses of optically thin H2 gas at 300 K and higher 
temperatures. We compare our results to a Chiang and Goldreich (1997, CG97) two-layer disk model of masses 0.02 M© and 
0.11 Mq. The upper limits to the disk's optically thin warm gas mass are smaller than the amount of warm gas in the interior 
layer of the disk, but they are much larger than the amount of molecular gas expected to be in the surface layer. If the two-layer 
approximation to the structure of the disk is correct, our non-detections are consistent with the low flux levels expected from 
the small amount of H2 gas in the surface layer. We present a calculation of the expected thermal H2 emission from optically 
thick disks, assuming a CG97 disk structure, a gas-to-dust ratio of 100 and Tgas = Tiust- We show that the expected H2 thermal 
emission fluxes from typical disks around Herbig Ae/Be stars are of the order of 10"'^ to 10"'^ erg s"' cm"^ for a distance of 
140 pc. This is much lower than the detection limits of our observations (5 x 10"'^ erg s"' cm"-). H2 emission levels are very 
sensitive to departures from the thermal coupling between the molecular gas and dust in the surface layer. Additional sources of 
heating of gas in the disk's surface layer could have a major impact on the expected H2 disk emission. Our results suggest that 
in the observed sources the molecular gas and dust in the surface layer have not significantly departed from thermal coupling 
(Teas/Tdust < 2) and that the gas-to-dust ratio in the surface layer is very likely lower than 1000. 
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1. Introduction 

Circumstellar disks surrounding low- and intermediate-mass 
stars in their pre-main sequence phase are the locations where 
planets presumably form. Such protoplanetary disks are com- 
posed of gas and dust. Their mass is initially dominated by gas 
(99%), specifically by molecular hydrogen (H2), which is the 
most abundant gas species. The dust constitutes only a minor 
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fraction of the total disk mass, however, it is the main source 
of opacity. Consequently, most of what we know observation- 
ally about protoplanetary disks has been inferred from studies 
of dust emission and scattering (for recent reviews see Henning 
et al. 2006; Natta et al. 2007; Dullemond et al. 2007; Watson et 
al. 2007). In order to understand the structure and evolution of 
protoplanetary disks, it is necessary to study their gaseous con- 
tent independently from the dust. For example, a basic physical 
quantity such as the disk mass is conventionally deduced from 
dust continuum emission at millimeter wavelengths assuming 
an interstellar gas-to-dust ratio of 100 (e.g., Beckwith et al. 
1990; Henning et al. 1994). If dust is physically processed in 
the disk, as should be the case in order to form planets, the gas- 
to-dust ratio must change with time. The disk dissipation time 
scale - another fundamental quantity required to disentangle 
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Table 1. Summary of the stellar properties and previous gas and dust observations. References: Acke et al. 2004 [A04], Acke 
et al. 2005 [AOS], Dent et al. 2005 [DOS], de Zeeuw et al. 1999 [D99], Dunkin et al. 1997 [DU97], Hipparcos catalogue [HIP], 
Pietu et al. 2003 [P03], van den Ancker et al. 1998 [V98], van der Plas et al. (in prep) [VP], Stempels & Gahm 2004 [S04]. 
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" Ages derived from the H-R diagram employing the pre-main-sequence evolutionary tracks of Palla and Stahler (1993). '' Total disk gas 
masses are derived from dust continuum emission at millimeter wavelengths assuming a gas-to-dust ratio of 100. ^ Group of each Herbig 
Ae/Be star disk according to the classification scheme of Meeus et al. (2001). Group I sources -flared disks- are sources with a rising mid-IR 
spectral energy distribution (SED) and Group II sources -self-shadowed disks- are sources with a flat SED (Dullemond et al. 2002). CTTS: 
Classical T Tauri star. Spectroscopic binary of period 2.4 days [S04]. ' We assume that both components have the same luminosity. ^ Disk 
mass by Jensen et al. (1996) scaled to our adopted distance of 145 pc. 



proposed giant planet formation scenarios (Pollack et al. 1996; 
Boss et al. 1998) - is deduced from observations of thermal in- 
frared excess emission produced by dust grains (Haisch et al. 
2001, 2005). Although recent studies (e.g., Sicilia-Aguilar et 
al. 2006) suggest a parallel evolution of the dusty and gaseous 
components, it still remains to be demonstrated that the gaseous 
disks disappear over the same time scale as the infrared excess. 

A variety of spectral diagnostics of the gas disk have been 
observed from the UV to the millimeter (see reviews by Najita 
et al. 2007; Dutrey et al. 2007). However, the only diagnostic 
that is potentially able to probe the warm gas in the regions 
where giant planets are thought to form is the mid-infrared 
(mid-IR) emission lines of H2. UV and near-infrared diagnos- 
tics only probe the innermost regions of the disk (R < few AU), 
and mm and sub-mm diagnostics are limited to probe the cold 
outermost regions of the disk (R >10 AU). 

Molecular hydrogen is by far the most abundant molecular 
species in protoplanetary disks. Unfortunately, H2 is one of the 
most challenging molecules to detect. Since H2 is a homonu- 
clear molecule, it lacks a permanent dipole moment and its 
transitions are thus electric quadrupole in nature. The small 
Einstein coefficients, characteristic of the quadrupole transi- 
tions, imply that H2 emission features are very weak. In ad- 
dition, in the case of protoplanetary disks, the H2 lines are 
not sensitive to the warm gas in the optically thick regions 
where the dust and gas are at equal temperature. Practical ob- 
servational challenges also have to be faced. The mid-IR H2 
emission from the disk needs to be detected on the top of a 
strong mid-IR continuum. From the ground, the mid-infrared 
windows are strongly affected by sky and instrument back- 
ground emission, and the H2 transitions at 12 and 17 jum lie 
close to atmospheric absorption lines highly dependent on at- 
mospheric conditions. The advent of high spectral resolution 
spectrographs mounted on larger aperture telescopes, allows 
for the first time the study of H2 emission from the ground, 
but the search is still limited to bright targets. From space, 
the problems of atmosphere absorption are alleviated and the 



7 = 2-0 feature at 28 fim is visible. However, the beam sizes 
are large and the spectral resolution of space mid-infrared spec- 
trographs are usually low when compared to ground-based fa- 
cilities, therefore, they are not very appropriate for small line- 
to-continuum ratios. 

H2 emission from protoplanetary disks in the mid-IR 
has been reported from ISO observations (Thi et al. 2001). 
However, subsequent ground-based efforts (Richter et al. 2002; 
Sheret et al. 2003; Sako et al. 2005) did not confirm the ISO de- 
tections. H2 emission in the mid-IR has been searched towards 
debris disks using Spitzer (Hollenbach et al. 2005, Pascucci et 
al. 2006, Chen et al. 2006) with no detection reported. Most 
recently, Bitner et al. (2007) and Martin-Zaidi et al. (2007) 
reported the detection of mid-IR H2 emission in two Herbig 
Ae/Be stars (AB Aur and HD 97048) from the ground, and 
Lahuis et al. (2007) reported the detection of mid-IR H2 emis- 
sion in 6 T Tauri stars with Spitzer 

Here, we report on a sensitive search for molecular hy- 
drogen emission from protoplanetary disks. We observed six 
southern nearby (d<400 pc) Herbig Ae/Be stars and one T 
Tauri star, employing the Very Large Telescope (VLT) imager 
and spectrometer for the mid-infrared (VISIR; Lagage et al. 
2004U ESO's VLT mid-infrared high-resolution spectrograph, 
and searched for H2 0-0 5(1) (7 = 3-1) emission at 17.035 jum 
and H2 0-0 5(2) (7 = 4-2) emission at 12.278 fim. The paper 
is organized as follows: in Sect. 2 we present the sample stud- 
ied, and the details of how the observations were performed. 
In Sect. 3 we discuss the data reduction. In Sect. 4 we deduce 
upper limits to the H2 fluxes, and using the optically thin ap- 
proximation, we derive upper limits to the mass of warm (150 
- 1000 K) gas in the disks. In Sect. 5 we discuss our results in 
the context of the Chiang and Goldreich (1997) two-layer disk 
model. Finally, we summarize our results and conclusions in 
Sect. 6. 



' http;//www.eso.org/instruments/visir 
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" The airmass interval is given from the beginning to the end of the observations. * Veiaj is the expected velocity shift of the spectra due to the 
reflex motion of the Earth around the Sun and the radial velocity of the star. ^ The standard stars were observed prior (P) and/or immediately 
following (F) the science observations. 



2. Observations 

Target selection. We selected a sample of well known nearby 
Herbig Ae/Be and T Tauri stars based on evidence of large disk 
reservoirs. The targets have either reported detections of cold 
CO gas at (sub)-mm wavelengths, or dust continuum emission 
at mm wavelengths. We chose stars with 12 jjm continuum 
fluxes > 0.5 Jy (otherwise too faint for acquisition with VISIR) 
and < 25 Jy (hard to detect weak lines on top of a strong con- 
tinuum). The physical properties of the targets are compiled in 
Table 1. 

Observations. The sources were observed in the first 
semester of 2006 and 2007 with VISIR, a combined imager and 
spectrograph designed for observations in the (a; 8 - 13 yum) 
and Q bands (=s 16.5 - 24.5 fim) (Lagage et al. 2004), mounted 
at the ESO-VLT Melipal telescope in Cerro Paranal, Chile. 
The H2 - 5(1) (7 = 3-1) line at 17.035 fim was ob- 
served in the high-resolution long-slit mode with a 0.4 arcsec 
slit, giving a spectral resolution R ^ 21000, or 14 km s The 
H2 - 5(2) (7 = 4 - 2) line at 12.278 fj.m was observed in 
the high-resolution echelle mode with a 0.4 arcsec slit, giving 
a spectral resolution R « 20000, or 15 km s~' . 

The total integration time in each line was Ih. The slit was 
oriented in the North-South direction. Sky background was 
subtracted by chopping the telescope by ~8" in the direction 
of the slit. Asymmetrical thermal background of the telescope 
was subtracted by nodding the telescope by ~8" in the direction 
of the slit. For correcting the spectra for telluric absorption and 
obtaining the absolute flux calibration, spectroscopic standard 
stars at close airmasses to that of the science targets were ob- 
served immediately preceding and following the 12 um expo- 
sure, and preceding or following the 17 jum exposurqj. Finally, 
we tried as much as possible to observe the stars at a moment 

^ The 12 fim window is much more affected by telluric absorption 
than the 17 fim window. In order to be able to select the standard star 
observation that provides the best telluric correction, observations of 
the H2 S(2) line require the observation of a telluric standard preceding 
and following the science observation. The 17 micron window is less 
aifected than the 12 micron window by telluric absorption, therefore. 



of high velocity shift between the telluric absorption feature at 
12.27 microns and the H2 S(2) line. This effect is due to the 
reflex motion of the Earth around the Sun. The typical PSF 
FWHM measured in the continuum was 0.5" (~4 pixels using 
the VISIR pixel scale of 0.127"/pixel) equivalent to a spatial 
resolution of 50 AU for a source located at 100 pc. A summary 
of the observations is presented in Table 2. 

3. Data reduction 

VISIR raw data consists of a collection of data cubes, in which 
each data cube is associated with data taken in one nod posi- 
tioiJE The odd planes in the data cube correspond to a half- 
chop cycle frame, the even planes correspond to the average 
of the difference of all preceding half-chop cycles (i.e. A,=Ai- 
Bi + A2-B2-i-...-i-A,-B,/i). The last plane corresponds to the av- 
erage of all half-cycle difference images. 

The first step in the data reduction process was to recover 
the individual half chop cycle frames from the data cube, and 
to guarantee that all of them have the same wavelength in the 
same row of pixels (in each 2D spectrum, the spatial direction 
is along the columns and the dispersion direction is along the 
rows). In each 2D half-chop cycle frame, at the position of the 
target spectrum, a cut with a size of 20 pixels in the spatial 
direction was made. From this cut, a spectrum of the atmo- 
spheric emission was obtained by the sum of the counts in the 
spatial direction (see uppermost right panel in Figs. 1 and 2). 
This atmospheric emission spectrum was rebinned by a factor 
of 50 (i.e. 50 pixels where previously we had 1). Using as ref- 
erence the atmospheric emission spectrum from the first half- 
chop cycle of the first nod, the ofi^set that each half-chop cycle 
raw frame required for having the center of the sky emission 
lines at the same pixel position as the reference frame was de- 
termined. The offset was given by the difference between the 



only one observation of a telluric standard (following or preceding the 
science observation) is required. 

^ The sequence of the observations is ABBA, where A and B corre- 
spond to two different nod positions. 
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Fig.l. Spectra obtained for the H2 0-0 S(2) (J=4-2) line at 12.278 yum. The left panel shows a zoom to the -100 to 100 km 
s ' interval of the atmospheric corrected spectra. A Gaussian of FWHM - 15 km s"^ and integrated line flux equal to the line- 
flux upper Umits obtained is overplotted at the expected velocity shifted location (vertical dashed Une, see Table 2). The central 
panel shows the full corrected spectra. Dotted lines show spectral regions strongly affected by telluric or standard star absorption 
features. The right panel shows the continuum normalized spectra of the standard star and the target before telluric correction. 
The uppermost right panel displays the sky spectrum from a half-chop cycle. The spectra are not corrected for the radial velocity 
of the targets. 



center of the Gaussian fitted to one sky emission line in the ref- 
erence spectra, and the center of the Gaussian fitted to the same 
sky line in the atmospheric emission spectrum of the half-chop 
cycle frame to be corrected. The typical offsets were a frac- 



tion of a pixel (before rebinning), and did not depend on the 
reference sky Une. Then, each 2D raw half-chop cycle frame 
was rebinned by a factor of 50 (i.e. 50 pixels where previously 
we had 1) in the dispersion direction, and it was shifted -in the 
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Fig. 2. Spectra obtained for the H2 0-0 S(l) (J=3-l) line at 17.035 yum. The left panel shows a zoom to the -100 to 100 km 
s~^ interval of the atmospheric corrected spectra. A Gaussian of FWHM - 15 km s"' and integrated line flux equal to the line- 
flux upper hmits obtained is overplotted at the expected velocity shifted location (vertical dashed hne, see Table 2). The central 
panel shows the full corrected spectra. Dotted lines show spectral regions strongly affected by telluric or standard star absorption 
features. The right panel shows the continuum normalized spectra of the standard star and the target before telluric correction. 
The uppermost right panel displays the sky spectrum from a half-chop cycle. The spectra are not corrected for the radial velocity 
of the targets. 
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dispersion direction - by the offset found. Finally, the obtained 
2D frame was rebinned to the original size of 256x256 pixels. 
This procedure corrected small differences in 2D half-chop cy- 
cles introduced by variations in the grating position within the 
instrument. 

Thereafter, with the corrected half-chop frames, the VISIR 
data cubes were re-built and processed with the VISIR pipeline 
(Lundin 2006). Bad pixels were detected in the half-cycle 
frames when their number of counts exceeded the limit of 
65000. They were cleaned by interpolation with neighboring 
pixels. In the case of the 17 jjm data (i.e. taken in the long-slit 
mode), the 2D frames were corrected for distortion by inter- 
polating the analytical optical distortion (see description in the 
VISIR user's manual, Siebenmorgen et al. 2006) and the value 
from the source pixels. The 2D frames in each nod position 
were created by averaging all the half-cycle difference images 
in the data cube. The nodding images were produced by av- 
eraging the image in the two nod positions A and (-B), and 
dividing the result by 2xDIT, where DIT is the detector inte- 
gration tim43- The nodded images are jittered and were thus 
shifted and added to form the final combined image using the 
offsets stored in the file header. 

The wavelength calibration was performed with the follow- 
ing procedure. The first half-cycle from the first nod position 
was "collapsed", producing a 1 -dimensional spectrum of the 
atmosphere. The centers of the atmospheric emission lines (in 
pixels) were found by fitting Gaussians to their profiles in a 
HITRAN model spectrum of Paranal's atmosphere emission 
(for details of the model, see Lundin 2006). The dispersion re- 
lation was found by a second degree polynomial fit between 
the centers of the atmospheric lines in the observed and model 
spectrum. 

The spectrum at each wavelength was extracted from the 
2D combined image by summing the pixels in the spatial direc- 
tion inside the PSF. In the case of the long-slit 17 fim data, one 
spectrum was extracted for each of the nod positions present in 
the final 2D combined image. The final extracted spectrum is 
the sum of these three spectra. The telluric correction was per- 
formed by dividing the extracted spectra by the extracted spec- 
tra of the standard star. Small shifts (of a fraction of a pixel) on 
the spectra of the standard star in the wavelength direction were 
applied until the best signal to noise spectra (i.e. best telluric 
correction) were found. The spectral flux calibration was made 
by multiplying the telluric-corrected science spectrum with a 
model spectrum (Cohen et al. 1999) of the standard stailfl 

In the case of the observations of the H2 S(2) line in HD 
34282, HD 100453 and HD 101412 the telluric standard stars 
showed photospheric absorption lines close to the location of 
the H2 line. In the case of the HD 34282 observation, the 
photospheric absorption is located outside the region where 
the H2 line is expected (after correction by the velocity shift). 
However, in the cases of HD 100453 and HD 101412 the ex- 
pected location of the S(2) H2 line partially overlaps with the 

The factor 2 is due to the fact that the on-source and the off-source 
images both contribute to the whole DIT. 

^ The catalogue of spectroscopic standard stars in the mid-IR is 
available at www.eso.org/VISIR/catalogue. 
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Fig. 3. Correction for photospheric absorption features in the 
spectrum of the telluric standard (STD) of HD 100453 and HD 
101412. In grey the observed STD spectrum, in black the cor- 
rected STD spectrum. The photospheric feature was removed 
by a Gaussian fit to the telluric line at 12.277 fim. The vertical 
dotted line displays the expected location of the S(2) H2 line in 
the science spectrum. 

location of the photospheric absorption feature in the spectrum 
of the standard star (see Fig. 3). This absorption feature is ob- 
served superimposed on the broad telluric absorption line at 
12.277 micron. To avoid the error induced by this feature in the 
telluric corrected spectra, we removed this absorption feature 
from the standard star spectra by a Gaussian fit to the underly- 
ing broad atmospheric absorption feature at 12.277 micron (see 
Fig. 3). After this correction, the science spectra were divided 
by the telluric standard and flux calibrated in the conventional 
way. 

4. Results 

In Figs. 1 and 2 we present a compilation of the obtained spec- 
tra. In the left panel, the atmosphere-corrected flux calibrated 
spectra in wavelength scale measured by the velocities in km 
s ' are shown. In the central panel, the full corrected spectra 
are presented. In the right panel, for reference, the normalized 
spectra of the science target and the standard star are displayed. 
None of the observed sources show evidence for H2 emission 
at 12 or 17 fim. 

Line flux upper limits. In the flux-calibrated science spectra, 
the standard deviation (cr) of the continuum flux was calculated 
in regions less influenced by telluric absorption, photospheric 
features of the standard star and close to the features of interest. 
In Figs. 1 and 2, the horizontal dotted line shows the 3cr limits 
of the derived continuum flux. 3cr upper limits to the integrated 
line flux were calculated by multiplying the 3o" flux with the 
instrument resolution Une widtlH of 15 km s In Table 3, a 
summary of the upper flux limits is presented. In Figs. 1 and 2 
a Gaussian line of FWHM -15 km s ' and an integrated flux 
equal to the upper flux limit is overplotted on each observed 
spectrum. The typical sensitivity limit of our observations is a 

^ Here we assume a pole-on disk and that the true-line width is nar- 
rower than the instrumental line width of 15 1cm s"' (i.e. the spectral 
line is not resolved). If the line is broader, the sensitivity limits in 
Table 3 decrease proportionally with the FWHM (in km s"') follow- 
ing F = F,,^^,-,* (FWHM/ 15). 



Table 3. Upper flux limits. 
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[Jy] 
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ergs s ' cm 2] [ log L/Lq] 
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[ log L/Lq] 
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a non-detection of 


the H2 S(l) and H2 S(2) lines with Icr upper limits of 0.9 x 10"'" erg s"' cm" 


^ assuming a line width of 500 km s ' . 
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Upper mass limits assuming optically thin H2 
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A„/ is the Einstein coefficient of the J — u - I transitiorQ and 
x„ is the population of the level u at the excitation temperature 
T in LTE; / is the conversion factor required for deriving the 
total gas mass from the H2-ortho or H2-para mass determined. 
Since Mh^ = MH^ortho + Mn.p^m, /=l+ortho/para for the S(2) 
line (a H2-para transition) and /=l-Hl/(ortho/para) for the S(l) 
line (a H2 -ortho transition). The equilibrium ortho-para ratio at 
the temperature T was computed using the Eq. 1 of Takahashi 
(2001). In Table 4, we present our results. 

The H2 S(l) and S(2) line fluxes constrain the amount of 
optically thin hot and warm H2 gas. The disks contain less than 
a few tenths of Jupiter mass of optically thin H2 at 150 K, and 
less than a few Earth masses of optically thin H2 at 300 K and 
higher temperatures. In Fig. 4 we present the mass of optically 
thin H2 as a function of the temperature derived from a S(l) 
and a S(2) line with a flux 10-'"* erg s-' cm-^ for a source 
at a distance of 140 pc. This mass is the minimum amount of 
optically thin H2 that would produce a detectable H2 line in our 
observations at a given temperatur^ We see that the S(l) line 

4.76 X 10-'" S-' for the S(l) line / = 3 - 1 transition, 2.76 x 
10-' S-' for the S(2) line 7 = 4-2 transition (Wolniewicz, Simbotin, 
& Dalgamo 1998). 

* We calculated a set of CG97 disk models around a prototypical 
Herbig Ae star (see Sect. 5, and parameters of MWC 480 in Table 5) 



line flux of 10-'** ergs s"' cm"^. Our flux limits are of the or- 
der of magnitude of the H2 S(l) and S(2) lines fluxes (1.1 and 
0.53 xlO-'** ergs s-' cm-^, respectively) reported by Bitner et 
al. (2007) for AB Aur and the H2 line fluxes of 0.33 - 1.70 
xlO-'"* ergs S-' cm-^ reported for the H2 S(2) line in the obser- 
vations by Lahuis et al. (2007). Here, we should note that the 
observations of AB Aur were performed with TEXES, which 
provides a spectral resolution of 100000, increasing the line-to- 
continuum contrast compared with our observations. 

Upper limits on the optically thin warm gas disk mass. 
Under the assumption that the H2 emission and the emission of 
the accompanying dust are optically thin, that the emitting H2 
is in local thermodynamical equilibrium, and that the source 
size is equal or smaller than VISIR's beam size, we derived 
upper limits to the H2 mass as a function of the temperature 
employing (Thi et al. 2001) 



Mgas =/ X 1.76 X 10-20 ^^'^ Mo, (1) 



where Fui is the upper limit to the integrated line flux, d is the 
distance in pc to the star, £„/ is the energy of the transition. 
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Fig. 4. Mass limits of optically thin H2 derived from H2 S(l) (solid line) and H2 S(2) (dotted line) as a function of the temperature 
for a line flux limit of 10"'"* ergs s"' cm"^ for a source at a distance of 140 pc. Dashed and dot-dashed lines show the gas mass as 
function of the temperature for a Chiang and Goldreich (1997) optically thick two-layer model for a low-mass (Mqisk = 0.02Mo) 
and a high-mass (Mdisk = 0.1 IMq) disk assuming a gas-to-dust ratio of 100. Mint is the mass of the interior layer. Msmf is the 
mass of the surface layer. H2 emission arises only from the optically thin molecular gas in surface layer of the disk. 



sets constraints for the mass of optically thin H2 at temperatures 
150 K < r < 500 K. For higher temperatures flie S(2) is a better 
trace:0. 

5. Discussion 

5.1. Comparison with the gas mass of a two-layer disk 
model 

To understand the mass upper limits derived, we suppose that 
the disk has a two-layer structure (Chiang & Goldreich, 1997). 
In this disk model, the infrared excess in the SED is produced 
by the radiation of small dust grains in the superheated opti- 
cally thin surface layer of the disk, in addition to the emission 
of cooler dust in the mid-plane. In the interior layer, the gas 
is thermalized by collisions with the dust grains, the column 
densities of H2 are larger than 10^^ cm"^ and the optical depth 
is higher than one (see Fig. 6). The interior layer is optically 
thick in both gas and dust even at mid-infrared wavelengths. 
Emission lines from an optically thick medium are only ob- 
servable when there are temperature gradients. The hot upper 
layer can provide emission lines against the cold mid-plane. 

Note that the two layer model is an approximation of the 
real structure of the disk. In fact, in the uppermost zone of the 
surface layer, the gas is thermally decoupled from the dust, and 
due to intense X-ray or UV heating may well be ionized and 
at temperatures of the order of 10'* K. This uppermost surface 

and found that the interior layer alR > 5AU at 17 and 12 yum would be 
optically thin if the mass of a disk with a size of 100 AU was smaller 
than 0.01 Mj, if a gas-to-dust ratio of 100 and Tg^s = Tdust is assumed. 

' For temperatures below 150 K the S(0) (7 = 2-0) line at 28 
^m provides better constraints to the optically thin gas than the S(l) 
line. However, as the Earth's atmosphere is completely opaque at these 
wavelengths, the S(0) line is only observable from space. 



is above the region where the hydrogen is molecular The H2 
emission arises from a layer below this uppermost surface at 
about Av=l, or a hydrogen column about 10^' to lO^-' cm"^ 
down from the surface, depending on the dust opacity. For the 
remainder of the discussion, one must keep in mind that when 
we speak about the disk surface layer, we refer specifically to 
the molecular gas layer of the disk's surface layer. 

Using a two-layer disk model implementation (CGplus, 
Dullemond et al. 2001) with physical parameters aimed to fit 
the SED of the prototypical Herbig Ae/Be stars, we computed 
the expected amount of gas in the interior and surface layer as 
function of the temperature for disks of mass 0.11 M© and 0.02 
Mq (see Table 5 and Fig. 50 In Fig- 4 we present our results. 
The mass limits derived from the H2 S(l) and S(2) line obser- 
vations are smaller than the amount of warm gas in the interior 
layer, but much larger than the amount of molecular gas in the 
surface layer. Fig. 4 shows that the amount of gas in the surface 
layer is very small (< 10"^ Mj ~ 3M®) and almost indepen- 
dent of the total mass of the disk. If the two-layer model is an 
adequate representation of the structure of the disk, the thermal 
flux levels 0///2 mid-infrared emission are below the detection 
limit of the observations, because the mass 0///2 in the surface 
layer is very small. 



In our models no puffed-up inner rim was used. We employed 
a model with a truncation radius at T = 5000 K. We assumed that 
Tgo.i=Trf„s, and used a gas-to-dust ratio of 100. For disk-models with 
TgojT^Tjus, see for example, Kamp and Dullemond (2004), Nomura 
and Millar (2005) and Jonkheid et al. (2006). 
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Fig. 5. Physical parameters of the two-layer disk model as a function of radius (R). Upper panels: high-mass disk model (Mdisk - 
O.IIM0, L»= 24.6L0). Lower panels: low-mass disk model (Mdisk - O.O2M0, L*= 2O.2L0). The indices int and surf refer to the 
interior and the surface layer. S is the surface mass density, T the temperature, M the disk mass and H/R the scale height. See 
Table 5 for details of the model parameters. 



Table 6. Observational efforts to detect mid-IR H2 emission from gas-rich protoplanetary disks. References: Bitner et al. 2007 
[B07], Cai-mona et al. (this work) [C07], Lahuis et al. 2007 [L07], Martin-Zaidi et al. 2007 [MZ07], Richter et al. 2002 [R02], 
Sako et al. 2005 [805], Sheret et al. 2003 [SR03], Thi et al. 2001 [TOl]. 



Instrument 


A/AA at 17yum 


LW 

[km s-i] 


beam 
[arcsec x arcsec] 


flux & flux limits 

[x 10-'* erg S-' cm-^] 


ref 


VISIR 


21000 


14 


0.4 X 0.4 


<0.5 


C07 




10000_£ 


30 


0.4 X 0.4 


2.4 


MZ07 


COMICS 


5000 ~ 


60 


0.6 X 0.4 


<1.0^ 


S05 


TEXES 


60000 


5 


2x2 


<1.4~ 


R02 




100000 


3 


0.4 X 0.4 




B07 


MICHELLE 


15200 


40 


1 X 1 


<6.73 


SR03 


ISO-SWS 


2000 


210 


14x 27 


<0.8 - 40 


TOl 


spitzer£ 


600 


500 


5x5 


0.32_f. 


L07 



" Observations of HD 97048 performed with a 0.75" slit, therefore, the different resolution. 

* The detection of Hj S(l) in HD 163296, MWC863, CQ Tau and LkCal5 by Thi et al. 2001 was not confirmed by Sako et al. 2005. 
The detection of H2 S(l) in CQ Tau and AB Aur by Thi et al. 2001 was not confirmed by Richter et al. 2002. 

Bitner et al. 2007 reported the detection of the H2 S(l), S(2) and S(4) lines in AB Aur. 
' The detection of H2 S(l) in HD 163296, AB Aur and GO Tau by Thi et al. 2001 was not confirmed by Sheret et al. 2003. 
^ From 76 sources studied one exhibits the H2 S(l) line and six the H2 S(2) line. None of the 8 HAeBes observed exhibit H2 mid-IR emission. 

* Detection of the H2 S(l) line in the T Tauri star Szl02. 



5.2. Derivation of expected S(1) and S(2) H2 emission 
from tlie two-layer disk model. 

Our observations are one of the most sensitive efforts to date 
to search for H2 emission from gas-rich protoplanetary disks 
(see Table 6). If the H2 emission is only produced from the 
upper layer of the disk, as the two layer model suggests, we 
would like to find out what flux levels we could expect from 
typical Herbig Ae/Be disks, assuming Tgas=T(iust and a gas-to- 



dust ratio equal to 100. For simplicity we first assume that the 
disk is viewed pole-on (/ = 0°). Here, we should note that the 
gas temperature can be higher than the dust temperature in the 
surface layer (e.g., Kamp and DuUemond 2004, Nomura and 
Millar 2005, Jonkheid et al. 2006, Glassgold 2007). 

In the two-layer disk model, at each disk radius the H2 
emission from the surface layer has to be observed on the top 
of the optically thick continuum of the interior layer, and the 
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Fig. 6. Column density of H2 and optical depth of H2 and dust as a function of radius R for the interior (upper panels) and the 
surface layer (lower panels) in the two-layer low-mass disk model (Mdisk - O.O2M0, Lt- 2Q.2Lq). A factor ten higher values 
in the H2 column density and optical depth of the interior layer are reached in the high-mass (Mdisk - O.IIM0, L^- 24.6Lo) 
model. See Table 5 for details of the model parameters. 



Table 5. Physical parameters of prototypical Herbig Ae/Be 
stars and their disks after Chiang et al. (2001). 



Parameter 


MWC 480 
high-mass 


HD 36112 
low-mass 


T. [K] 


8890 


8465 




2.1 


2.1 


M, [Mo] 


2.3 


2.2 


L. [he] 


24.6 


20.2 


d [pc] 


140 


150 


Eo [g cm"-] 


8000 


1000 


ao [AU] 


100 


250 


H/h 


1.7 


1.5 


P 


1.5 


1.5 


Mdisk [Mq] 


0.11 


0.02 



optically thin continuum of the dust at the surface layer. The 
total intensity (line plus continuum) emitted at each radius (R) 
is given by 

Iy(R) = Bv.(ri„,)e-"'- + Bv(rsurf)(l - e"^-»0. (2) 

Here, By(7'int)e"'^'™'f is the contribution of the interior layer to 
the intensity. By(Ti„i) is the black body emission of the interior 
layer at temperature Tint seen through the surface layer optical 
depth (Tv,„,j = Tv,„,„„,j -1- Ty^,,,,J . The contribution of the surface 



layer is Bv(TsuYt)(i - e~'^'<"rf). It is the optically thin emission of 
the gas and dust in the surface layer at temperature Tsmi- 
The dust optical depth in the surface layer is given by 

''"vdiislsLirf ~ '<'vdusi^dusts„rf ; (3) 

where Ky^^^^ is the dust mass absorption coefficient at the fre- 
quency V, and £dusts„rf is the dust surface mass density in the 
surface layer of the disk given by the CG97 model. 

The H2 gas optical depth in the surface layer at the frequen- 
cies V in the vicinity of the H2 line at frequency vq is 

T"vga.„„,r = fvii„2:ga.s,„,-/2, (4) 

with 2ga,s,„^| being the gas surface mass density in the disk sur- 
face layer (equal to 100 times 2dust,„,fl3. and /2 the conversion 
factor from the total gas surface density Sgass^rr to the H2-ortho 
or H2-para surface densities. Since Mgas - Monho + Mp^^d, then 
./2=l/(l+ortho/para)for the S(2) line (a H2-para transition) and 
/2=l/(l-i-l/(ortho/para)) for the S(l) line (a H2-ortho transi- 
tion). The gas mass absorption coefficient (/Cv,;^) in the vicinity 
of the H2 transition at frequency vo is given by 

'^vii„, = i^'^ixiBu, - x„Bui)(py. (5) 



" We used the opacities by Laor & Draine (1993). 
Here we assume that all the gas is H2. 
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Fig. 7. Surface optical depth and predicted H2 S(l) and S(2) emission line flux for the regions where the surface temperature (Ts) 
~ 150K (upper panel) and ~ 300K (lower panel) in the two-layer high-mass disk model (Mdisk = O.IIM0, L,= 24.6Lq). T, is the 
interior temperature. An instrumental line FWHM of 15 km s"' was assumed. Note that the flux scale in the line flux-continuum 
plot is 10 times smaller in the S(2) line plot as in the S(l) line plot. 



Here, Inip is the mass of the H2 molecule (2 times the mass of 
the proton nip), xi and jc,, are the lower and upper level popula- 
tion of the transition J - u - I, and B;„ and B„; are the proba- 
bilities of absorption and emission. 

We assumed that the gas is in LTE and derived the level 
populations using the Boltzmann equation. The line profile (py 
is assumed Gaussian with a FWHW^\ equal to 15 km s"' : 

1 I '-'-0 \- 

(pv - ]=e ^ '^y ' , (6) 

cTv yjn 

where o-y = (vo/c) * FWHM I {l \^). The expected flux 
dFy(R) at each radius (R) in the disk is given by 

2nRdR 

dFy(R) = iyiR)dn = 

with dQ being the solid angle, and dR the width of the region 
at a distance R from the central object where the surface tem- 
perature is equal to Tsurt and the interior temperature is equal 
to Tint. The distance to the source is d. 

Fig. 7 shows two examples of the surface optical depth 
(Tv^,,^i.) and the continuum subtracted S(l) and S(2) H2 predicted 
line fluxes for the regions with surface temperature ~150 K 
(R -41 AU upper panel) and -300 K (R ~7 AU lower panel). 
We assumed for the calculation a face-on high-mass disk model 

We ciieclced tliat the optical deptli is smaller than one with a (p 
with thermal width cr,/, = ^kTJ{2mp), therefore, could be assumed 
from the start with a width equal to the instrument resolution. 



at a distance of 140 pc. The column density of H2 in the inte- 
rior layer of the disk is larger than 10^^ cm"^. The total optical 
depth peaks at the H2 S(l) and S(2) rest wavelengths. 

The contribution to the total flux in each line varies as a 
function of the distance to the central source. Since at each ra- 
dius, the surface layer optical depth of dust and H2 gas is much 
smaller than one (see Fig. 6) Eq. 2 can be approximated by 

Iy{R) = By{Ti,d + Bv(rsurf)rv,„,„„,r + ■Bv(rsu,-f )Tv,„,„r • (8) 

In this way, the contributions to ly by the interior layer Bv(7'int), 
by the dust in the surface By(Tsud) ■'"vj„s,.„,f and by the H2 line 
Bv(7^surf) Tvju^.irr at the surface can be separated at each radius. 
Combining Eq. 7 with Eq. 8, we can disentangle the different 
contributions to the flux dFy as a function of the radius 

dFy(R) = dFyJR) + c/Fv,„„.„,r(^) + dFy^.,^^JR). (9) 

In the left panels of Fig. 8 we show the contributions, 
in percentage of the peak flux at each radius, by the inte- 
rior layer {dFy.^XR)/dFy{R), solid line), the dust in the sur- 
face (dFy^^^^^^^^^^(R)/dFy(R), dashed lines) and the H2 emission 
(dFy ^^^^^^^(R) /dFy{R), dotted lines). Up to a few AU, the main 
source of continuum emission is the dust in the interior layer; 
at larger radii, the principal source of continuum emission is 
the superheated dust in the surface layer. The H2 line contri- 
bution is always very small, typically less than 1% of the total 
continuum flux. 

In the central panels of Fig. 8 we show the evolution of line 
peak flux minus the continuum as a function of radius. In the 
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Fig. 8. Total S(l) and S(2) line fluxes for the two-layer high-mass disk model (upper six panels) and low-mass disk model (lower 
six panels). Left panels present the percentage contribution to the peak flux by the optically thick interior layer (solid line), the 
optically thin dust in the surface layer (dashed line) and the H2 line (dotted Une) as a function of radius for a pole-on disk {i=0 °) . 
Central panels show the peak flux minus the continuum as a function of radius for pole-on disk. Right panels display the expected 
total Ho line flux for three disk inclinations: nole-on fsolid lineV i=^0 ° fdashed line^ and i=60 ° fdotted lineV 
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case of the S(l) line the maximum contribution to the total flux 
is given by the region around 50-1 00 AU. In the case of the S(2) 
line the most important contribution to the line flux is made by 
material around 5-20 AU. 

The total expected flux from the disk is the sum of the con- 
tributions of each radius. In the right panels of Fig. 8 we present 
the total expected line from the disk. The line strength over the 
continuum is less than 0.05 Jy for the S(l) line and less than 
0.005 Jy for the S(2) line. The total integrated flux of the line is 
of the order of lO""" erg s"' cm"^ for the S(l) line and 10"'^ erg 
s"' cm"^ for the S(2) line. These line flux levels are two orders 
of magnitude below the sensitivity limits of our observations. 

We note that our H2 line flux predictions from the two-layer 
model only take into account thermal excitation of the gas, 
therefore, show the minimum amount of expected flux. If ad- 
ditional excitation mechanisms, such as X-rays or UV heating, 
are taken into account, the line fluxes are expected to be higher 
because the gas is hotter In sophisticated models of the disk 
around the T Tauri star TW Hya (d = 56 pc) by Nomura et al. 
(2005, 2007), the predicted H2 S(l) and S(2) line fluxes are of 
the order of a few 10"'^ erg s"' cm"^ and the line-to continuum 
flux ratios are approximately 10"^. For a star at 140 pc these 
line fluxes translate to a few 10"'^ erg s~' cm"^, fluxes simi- 
lar to our two-layer modeQ These theoretical H2 line fluxes 
are still two orders of magnitude below the detection limits of 
our observations. In addition, the achievement of line-to contin- 
uum flux ratios < lO^^* is a challenging task for present ground- 
based observations. In our observations, for example, the best 
Une/continuum ratio achieved is 10 

5.2.1 . The effect of the disk inclination 

Until now, we assumed for simplicity that the disk is viewed 
pole-on. Now, if the disk is inclined by an angle /, the line pro- 
file <pv used in Eq. 6 will be doppler-shifted for each parcel of 
gas that has a velocity component in the line of sight. For a par- 
cel of gas located at radius R and at azimuth^ 9, orbiting a star 
of mass Mi,, the doppler shift Av due to the Keplerian motion 
of the gas is 



Av - (vo/c) sin(0)sin{i) y ^ * , (10) 

with c being the speed of light. The doppler-shifted line profile 
<pv shifted will be 

<^v shifted = V <-v (11) 

Using Eq. 1 1 instead of Eq. 6 and employing Eqs. 3 to 8, the in- 
tensity I(R, 6, i) and the flux dFyiR, 6, i) emitted by each parcel 
of fluid of solid angle (R dR d9)/d^ is calculated. The expected 
total emission by the disk is the sum of the contributions of all 
the fluid parcels from to R,„a.x, the inner and outer radius of 

Note that T Tauri stars are less luminous than Herbig Ae/Be stars. 
The fluxes are similar because the Nomura models include additional 
sources of heating beyond dust-gas collisions. 

We define 6=0 as the angle where the velocity component of the 
motion is in the line of sight. 



the emitting region. Assuming that /?„„■„ = 0.5 AU and that R^ax 
is the disk outer radius, we calculated the expected H2 S(l) and 
H2 S(2) line profiles for the high-mass and the low-mass disk 
at inclinations of 30 and 60 degrees. Our results are presented 
in the right panels of Fig. 8. The net effect of the inclination is 
the broadening of the line and a decrease of the peak flux. The 
S(2) line is more affected than the S(l) line by the inclination. 
This is due to the fact that S(2) comes from warmer gas, which 
is closer in and therefore orbiting faster. Since most of the ma- 
terial responsible for the S(l) and S(2) emission is located at 
distances larger than a few AU (see central panels of Fig. 8) 
even at high inclinations, at our spectral resolution, the double 
peaked profile is not exhibited. This result suggests that higher 
temperature lines emitted in regions closer to the star (e.g., the 
S(3) or S(4) lines) are more suitable for the investigation of 
inner disk dynamics. 

5.3. Detections and non-detections of H2 mid-IR 
emission from protoplanetary d/sks 

Our calculations of the expected thermal H2 line fluxes from 
CG97 disks around typical Herbig Ae/Be stars, and sophisti- 
cated models of H2 emission from disks around T Tauri stars 
by Nomura et al. (2005, 2007) show that the expected emission 
fluxes of H2 from optically thick disks in the mid-lR are very 
weak (< 10"'^ erg s"' cm"^ for a source at 140 pc). These mod- 
els provide a good explanation for the numerous non-detections 
of H2 mid-lR emission from optically thick disks reported so 
far (GG Tau, HD 163296, Richter et al. 2002; CQ Tau, Sheret et 
al. 2003; L551 1RS5, DG Tau, GW Ori, Richter et al. 2004; HD 
163296, MWC 863, CQ Tau and LkCa 15, Sako et al. 2005; 62 
T Tauri and 8 Herbig Ae/Be stars, Lahuis et al. 2007; UX Ori, 
HD 34282, HD 100453, HD 101412, HD 104237, HD 142666 
and HD 319139, this work). 

However, some detections of mid-IR H2 emission from 
disks have been reported in the Herbig Ae/Be stars AB Aur 
(Bitner et al. 2007i3 and HD 97048 (Martin-Zaidi et al. 
20070 and in the T Tauri stars Sz 102, EC 74, EC 82, Ced 
1 10 IRS6, EC 92, and ISO-Cha 237 (Lahuis et al. 20070 An 
interesting question to address is the reason for the high H2 
fluxes observed in these sources. Bitner et al. 2007, based on 
the gas temperature deduced from the observed lines (670 K) 
and constraints on the line emitting region (18 AU) by a line 
profile fit, concluded, by comparison to a disk model of AB 
Aur (Dullemond et al. 2001), that in order to explain the ob- 
served emission, an additional heating mechanism (X-rays or 
UV heating, e.g., Glassgold et al. 2007; Nomura et al. 2007) 
is needed to heat the gas in the disk upper layers (i.e. Tg^ssurf/ 
Tdustsuif > !)■ A second possibility to explain high H2 flux lev- 
els is to invoke a gas-to-dust ratio much larger than the canon- 
ical value of 100 in the surface layer of the disk (for example 
due to dust coagulation and sedimentation). Martin-Zaidi et al. 
(2007) estimated the gas-to-dust ratio in HD 97048 from the 

H2 S(l), S(2) and S(4) lines. 
" H2S(l)line. 

H2 S(2) line, in addition the H2 S(l) and S(3) lines in Sz 102 and 
the H2 S(3) Une in EC 74. 
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Fig. 9. H2 S(0), S(l), S(2), S(3), S(4) line fluxes for the two-layer low-mass disk model as a function of the gas-to-dust ratio in 
the surface layer for Tgas/Tdust in the surface ranging from 1.0 to 2.0. The dotted horizontal Une presents typical ground detection 
limits. The S(0) line is only observable from space. 



detected H2 S(l) line and the dust mass required for producing 
the observed flux level in the continuum. They obtained gas- 
to-dust ratios ranging from 3000 to 14000 for gas temperatures 
between 150 K to 1000 K. 

In order to explore the influence of a change in the gas-to- 
dust ratio and the thermal decoupling of gas and dust in the 
surface layer, we calculated the expected S(0) to S(4) H2 line 
fluxes from the low-mass disk model as a function of the sur- 
face gas-to-dust ratio for Tgassurf/Tdustsurf ranging from 1.0 to 
2.0. We proceeded to use a simplified prescription. The output 
of a CG-plus model was taken as the disk surface mass den- 
sity and temperature of the dust in the upper and interior layer. 
To obtain the molecular gas surface mass density in the upper 
layer, the dust surface mass density in this layer was multiplied 
by the gas-to-dust ratio. To obtain the molecular gas tempera- 
ture in the upper layer, the dust temperature in this layer was 
multiplied by Tgassurf/Tdustsurf • Our results are presented in Fig. 
9. 

As expected, an increase in the gas-to-dust ratio or 
Tgassurf/Tdustsurf yields an increase in the H2 flux emission lev- 
els. However, deviations from T'gassurf=7dustsurf have a much 
larger impact on the increase of the strength of the H2 lines 
than the changes in the gas-to-dust ratio. For example, to ob- 
tain flux levels over 10"^"* erg s"^ cm"^ in the S(l) line - suffi- 
cient to explain the detections by Bitner et al. and Mantin-Zaidi 
et al. - it is just required that Tgassurf/Tdustsurf > 2.0. In contrast, 
to achieve similar S(l) flux levels the gas-to-dust ratio should 
increase by a factor of thirty. The most Ukely scenario is that in 
the molecular upper layer of the disk, both the gas temperature 
and the gas-to-dust ratio change. For example, a simultaneous 
increase in the gas-to-dust ratio (> 1000) and gas temperature 
(Tgas surf /T dust surf > 2.0) can boost the flux levels of the H2 S(2) 
line to detectable levels. Therefore, the non detections in our 
observations suggest that in the observed disks, the molecular 



gas and dust in the surface layer has not significantly departed 
from thermal coupling and that the gas-to-dust ratio in the sur- 
face layer is not larger than 1000. 

One interesting aspect is that observations from space 
(Lahuis et al. 2007) report more numerous detections of the 
S(2) line than detections of the S(l) Une, although in our LTE 
models the S(l) line is in general stronger than the S(2) line. 
Note that for a given temperature the amount of ortho H2 (re- 
sponsible for the S(l) Une) is higher than the amount of para 
H2 (responsible of the S(2) line). This discrepancy suggests 
that the level populations and the ortho/para ratio of the H2 
producing the Une in the observed sources are not in LTE. One 
additional piece of information suggesting that the H2 gas re- 
sponsible for the Une emission is heated by a mechanism other 
than collisions with dust, is the detection of a H2 S(4) line 50% 
stronger than the H2 S(l) Une in AB Aur (Bitner et al. 2007). In 
our LTE disk models the H2 S(4) line is always much weaker 
than the H2 S(l) Une. Our results confirm the suggestion by 
Bitner et al. (2007) that in AB Aur the H2 emission is produced 
in a narrow region (at 18 AU), in which the gas has been ex- 
cited to a high temperature (T ^670K). We speculate that this 
could be the location of shock-heated gas in the AB Aur disk. 



6. Summary and conclusions 

We observed a sample of nearby pre-main sequence stars with 
evidence for cold {T < 50 K) gas disk reservoirs and searched 
for emission of the warm gas with T > 150 K, which is ex- 
pected to be present in the inner region of these disks. None of 
the targets show any evidence for H2 emission at 17.035 yum or 
at 12.278 yum. From the 3cr upper limits of the H2 line fluxes, 
we found stringent upper limits to the mass of optically thin 
warm H2. The disks contain less than a few tenths of Jupiter 
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mass of optically thin H2 at 150 K, and less than a few Earth 
masses of optically thin H2 at 300 K and higher temperatures. 

Assuming that Tgas=Tjiust and a gas-to-dust ratio of 100, we 
compared our results to models of disks employing a Chiang 
and Goldreich (1997) optically thick two-layer disk model of 
masses 0.02 M© and 0.11 M©. The upper Umits to the disk op- 
tically thin warm gas mass are smaller than the warm gas mass 
in the interior layer of the disk, but they are much larger than 
the amount of molecular gas expected to be in the surface layer. 
The amount of mass in the surface layer is very small (< 10"^ 
My) and almost independent of the total disk mass. We calcu- 
lated the expected H2 S(l) and H2 S(2) Une fluxes emitted from 
a two-layer disk for the low-mass and the high-mass cases as- 
suming a distance of 140 pc, and LTE thermal emission. The 
predicted hne fluxes of the two-layer disk model are of the or- 
der of ~10~'^ - 10"^^ erg s~' cm"^, much smaller than the de- 
tection limits of our observations (5 x 10"'^ erg s"' cm"^). 

If the two-layer approximation to the structure of the disk 
is correct, we are essentially "blind" to most of the warm H2 
in the disk because it is located in the optically thick interior 
layer of the disk. Our non-detections are explained because of 
the small flux levels expected from the little mass 0///2 present 
in the optically thin surface layer. Naturally, the two-layer disk 
model is only an approximation of the real structure of a pro- 
toplanetary disk. The pufl'ed-up inner rim, which could be an 
important contributor to the H2 emission, is not included in 
our models, and in reality there will be a smooth transition 
zone between the disk hot surface layer and the cool disk inte- 
rior, which again could contribute significantly to the H2 emis- 
sion. In addition, the surface layers of the disk are hkely to 
have gas temperatures hotter than the dust temperatures or a 
gas-to-dust ratio higher than the canonical value of 100. Both 
eff'ects could potentially increase the H2 emission by signifi- 
cant amounts. We presented additional calculations in which 
the gas-to-dust ratio and the temperature of the molecular gas 
in the surface layer of the disk were increased. We showed that 
detectable S(l) and S(2) H2 Une flux levels can be achieved if 
Tgas surf/Tdust surf > 2 and if the gas-to-dust ratio in the surface 
layer is greater than 1000. H2 emission levels are very sensitive 
to departures from the thermal coupling between the molecu- 
lar gas and dust in the surface layer Our results suggest that 
in the observed sources the molecular gas and the dust in the 
surface layer have not significantly departed from thermal cou- 
pling and that the gas-to-dust ratio in the surface layer is very 
likely lower than 1000. A definitive interpretation of our results 
awaits the development of future, more sophisticated models. 
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